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COMPUTER PROGRAM FOR PRELIMINARY DESIGN ANALYSIS 
OF AXIAL -FLOW TURBINES 
by Arthur J. Glassman 
Lewis Research Center 

SUMMARY 

This report presents a computer program for the preliminary design analysis of 
axial -flow turbines. The computations are based on mean-diameter flow properties and 
do not consider any radial gradients. Given as input to the prc^ram are power or pres- 
sure ratio, mass flow rate, inlet temperature and pressure, rotative speed, inlet and exit 
diameters (either hub, mean, or tip), exit radius ratio or stator exit angle, turbine loss 
coefficient, and gas properties. Computations are then performed for any specified 
number of stages and for any of three types of velocity diagrams (symmetrical, zero 
exit swirl, or impulse). Exit turnir^ vanes can be included in the design. The pregram 
output includes inlet and exit annulus dimensions, exit temperature and pressure, total 
and static efficiencies, blading angles, and last-stage critical velocity ratios. 

The analysis method, a complete description of input and output, and a FORTRAN IV 
program listing are presented in this report. Sample cases are included to illustrate 
use of the program. 


INTRODUCTION 

The preliminary analysis of a power or propulsion system involves many repetitive 
calculations to determine system performance, component performance, and component 
geometries over a range of conditions. This must be done in order to eventually deter- 
mine the best system and operating conditions. For this type of screening analysis, 
complete design accuracy and detail for the components are not necessary . Approximate 
and rapid generalized procedures rather than complex and time-consuming detailed de- 
sign procedures are sufficient to yield the desired component overall geometry and per- 
formance characteristics. 

This report presents a computer program for the preliminary design analysis of 


axial-flow turbines. The analysis is based on mean-diameter flow properties and does 
not consider any radial gradients. Input design requirements include power or pressure 
ratio, mass flow rate, inlet temperature and pressure, and rotative speed. The design 
variables include inlet and exit diameters, stator angle, and number of stages. Compu- 
tations are jjerformed for any of three types of velocity diagrams (symmetrical, im- 
pulse, or zero exit swirl) by assuming the same shape diagrams for each stage. The 
program output includes inlet and exit annulus dimensions, exit temperature and pres- 
sure, total and static efficiencies, blading angles, and last-stage velocity ratios. 

The analysis method, a complete description of input and output, and a FORTRAN IV 
pregram listing are presented in this report. Sample cases are included to illustrate 
use of the pregram. 


SYMBOLS 


A stage loss parameter 

2 2 

Aan annulus area, m ; ft 

B exit loss parameter 

C blade loss parameter 

dimensional constant, 2ii rad/rev; 60 sec/min 
Cg dimensional constant, 1; 550 ft-lb/(sec)(hp) 

Cp heat capacity , joules/(kg)(K); Btu/(lb)(^) 

D diameter, m; ft 

E squared ratio of stage-exit axial velocity to stage-average axial velocity 
F blade loss weighting factor 

g dimensional constant, 1; 32.2 ft-lbm/(sec^)(lbf) 

Ah specific work, joules/kg; Btu/lb 

i stage number i, i = l,2, . .'.,n 

J dimensional constant, 1; 778 ft-lb/Btu 

K turbine loss coefficient 

M Mach number 

N rotative speed, rad/sec; rpm 

n number of stages 

P shaft power, watts; hp 
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p pressure, N/m^; Ib/ft^ 

R gas constant, joule s /(kg) (K); ft-lbf/(lbm)(^) 

Re Reynolds number 

r radius, m; ft 


temperature, K; °R 
blade speed, m/sec; ft/sec 
absolute gas velocity, m/sec; ft/sec 
relative gas velocity, m/sec; ft/sec 
mass flow rate, kg/sec; Ib/sec 
absolute-flow angle from axial direction, deg 
relative -flow angle from axial direction, deg 
heat capacity ratio 
static efficiency 
total efficiency 
speed-work parameter 
viscosity, (N)(sec)/m^; lb/(sec)(ft) 
density, kg/m^; Ib/ft^ 

Subscripts; 
a first stage 

cr critical 

ev exit vane 

ex turbine exit 

h hub 

i stage i, i = 1,2, . . . , n 

in turbine inlet 

m mean section 

n last stage 

ro rotor 

st stator 

t tip 


T 

U 

V 
w 

w 

Of 

y 

V 

n' 

X 

4 

p 
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u tangential component 
X axial component 

1 stator exit 

2 rotor exit 
Superscripts: 

turbine overall 

' absolute total condition 

" relative total condition 


METHOD OF ANALYSIS 

The method is based upon an analysis of the flow at the turbine mean diameter. 
Radial gradients of the flow properties are not considered. Specific heat ratio is as- 
sumed constant throughout the turbine. For any given tvirbine, all stages, except the 
first, are specified to have the same shape velocity diagram. The first stage differs 
only in that the inlet flow is axial. The velocity diagram shape depends upon the speed- 
work parameter value and the specified type of velocity diagram. Three types of veloc- 
ity diagram are considered: symmetrical, zero exit swirl, and impulse. These three 
types of velocity diagram are shown in figure 1 for three values of speed-work param- 
eter. 

Various input options dictate the exact nature of the calculation procedure. There 
is, however, one basic procedure that is direct and without iteration. This basic pro- 
cedure will be presented and then the alternate procedures required for the various in- 
put options will be discussed. The computations can be done either in SI units or in U.S 
customary units. 


Basic Calculation Procedure 

The required inputs for the basic procedure are shaft power P, mass flow rate w 
inlet total temperature T:^, inlet total pressure p|^, rotative speed N, inlet mean 
diameter exit mean diameter stator exit angle a^, gas constant R, 

specific heat ratio y , viscosity p, loss coefficient K, and squared ratio of stage-exit 
to stage -average axial velocities E. Also specified for each calculation are the number 
of stages n and the type of velocity diagram. For a multistage turbine, the input var- 
iable specified as inlet diameter is used to calculate first rotor blade speed and annulus 
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dimensions at the first rotor exit. Therefore, it is truly an inlet diameter only if the 
hub and tip diameters are assumed constant across the first stage. For a one-stage 
turbine, the specified exit diameter is used for the calculations, and the inlet diameter 

is of no significance. 

The first- and last-stage blade speeds are 


ttND 


U 


m, in 


n r- 


( 1 ) 


( 2 ) 


For more than two stages, it is assumed that the stage blade speeds vary linearly be- 
tween the first- and last-stage values. Therefore, 

Un - Ua /ON 

U. = ^ (i - 1) + U„ (3) 

^ n - 1 


Turbine specific work is 



J w 


and is equal to the sum of the specific work of the stages: 

n 

Ah’ = ^ Ah! 
i=l 

Expressing stage specific work in terms of stage speed-work parameter as 


(4) 


(5) 


U" 

Ah! = — 

' gJX 

and substituting equation (6) into equation (5) yield 


( 6 ) 
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( 7 ) 


Ah' = 


n 



i=l 


Since the velocity diagram shape is specified to be the same for all stages, the speed - 

work parameter is the same for all stages and is computed by rearranging equation (7) 
as 

gJAh' 


The value of X is the primary factor determining timbine efficiency . 

The method used for computing turbine static efficiency is basically similar to that 
presented in reference 1, but has the following additional featvires: (1) the turbines con- 
sidered in this report are not restricted to a constant mean-section diameter, (2) exit 
vanes to provide axial flow leaving the turbine can be included in the design, and (3) the 
velocity diagrams can be specified to be symmetrical. The efficiency computation 
method is explained fully in reference 1, and only the key equations are presented in 
this section. 

With turbine reheat neglected, turbine static efficiency can be expressed as 


7 = 


Ah' 


n-1 


( 9 ) 


Ah!^ „ + 


'id,a " XI ^^id,i 
i=2 


+ Ah. , 
id,n 


Dividing numerator and denominator by Ah' and introducing stage efficiencies yield 


V = 



i=2 


(10) 


Dividing equation (6) by equation (7) shows that 
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( 11 ) 


Ah! 
1 

Ah' 


Uf 



i=l 


Substituting equation (11) into equation (10) and recc^nizing that 



finally yield 


1 


1 ^ +1- 


n 

+ 

i=l 

[ ¥ 

m 

i=l 


(12) 


(13) 


The stage-total and last-stage -static efficiencies are 





(14) 


and 




n 


X 


X + i (A + B) 

2 


(15) 


The stage loss parameter A is expressed as 
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(16) 


KRe 


- 0.2 


cot a 




ro ro 


ev' 


1 


The constant of proportionality K, called the turbine loss coefficient in this report, 
must be determined empirically. On the basis of comparisons of predicted with experi- 
mental efficiencies, a value of K = 0.4 was selected in reference 1. For large tur- 
bines of recent airbreathing engines, a value of K = 0.35 seems better. The Reynolds 
number used in this calculation is defined as 


Re 


2w 


pD 


m,in 


(17) 


Some of the terms within the parentheses in equation (16) are the same for all 
cases, while others depend on stage location, velocity diagram type, and use of exit 
vanes. The rotor weighting factor and rotor loss parameter are the same 

for all cases: 


F =2 
ro 


(18) 




(19) 


For all stages other than last stages and for last stages where exit vanes are not used, 
the exit vane loss parameter is 


C =0 
ev 


( 20 ) 


For last stages of turbines having exit vanes. 


AVJ I 


(21) 


Axial inlet flow is assumed for all first-stage stators, for which the stator loss param- 
eter is expressed 
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2 


( 22 ) 


C 


St 


1 + 


2 cot^ff- 



For all stators other than first-stage stators, 


C 


st 



+ 



(23) 


The stator weighting factor also depends on whether or not the stator is a 

first-stage stator and further depends on the type of velocity diagram. The inlet and 
exit swirl parameters ^/AV^ and equations (19) to (23) also depend 

on the type of velocity diagram. The following table presents the relations for evaluating 
the stator weighting factor and the swirl parameters: 


Stage 

Velocity 

Stator 

Inlet swirl 

Exit swirl 


diagram 

weighting 

parameter. 

parameter, 


type 

factor, 


V o/AV 
u, 2 u 



fst 



First 

Symmetrical 

1 

X + 1 

X - 1 

2 

2 


Zero exit swirl 

1 

1 

0 


Impulse with 

1 

X.i 

x-i 


X 0 . 5 


2 

2 


Impulse with 

1 

x.i 

x-i 


X 0.5 


2 

2 

Intermediate 

Symmetrical 


X + 1 

X - 1 

2 — X 

2 

2 

and last 

Zero exit swirl 

1 

1 

0 


Impulse with 

2(1 - X) 

X +i 

x-i 


X < 0.5 


2 

2 


Impulse with 

1 

X +- 

x-i 


X ^ 0. 5 

1 

2 

2 


9 




The exit loss parameter B of equation (15) is expressed as 


B = E cot^o 




u,l 


1 


AV 


U/ 



where 


=!h:2 

I AV / AV 

\ "4v u 


with no exit vanes and 



= 0 


with exit vanes. 

The turbine exit velocities and state conditions are computed as follows: 


V 

u,ex 


U 

AV = — 
u,n ^ 


u,l,n 


- u,l 
AV, 


AV, 


u.n 


V „ = AV 

u,2,n u,n 

u 


^x,n = 


V = J~E V 
x,ex T ''x,n 


V 


,1 9 « no exit vanes 

u, /S, n 


( 24 ) 


(25) 


(26) 


(27) 

(28) 

(29) 

(30) 

(31) 

(32) 
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0 


exit vanes 



V = V 

''ex r ''x,ex *u,ex 


y R 


P V - 1 J 


/ — \r/(r-i) 

p ,p: 


Ah’ 

T' = T’ - ’’ 

ex in _ 

P 


T = T' 
ex ex 


Pex ^ex I -p 


, \y/(y-i) 


In order that the turbine total and static efficiencies be consistent with the computed 
exit velocity, the total efficiency is computed as 


c T! 1 
p in 


, \(y-l)/V 


Exit annulus area, radius ratio, and hub and tip diameters are obtained as follows: 


an, ex y 


ex x,ex 
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an, ex 



1 - 


ttD' 


m,ex 


1 + an, ex 


ttD" 


m,ex 


D 


t,ex 



D 


h,ex 



( 42 ) 


( 43 ) 


( 44 ) 


Absolute and relative flow angles, which are the same for each stage, are computed 
from the last-stage velocities: 


W 

W 


1 V „ 

. tan-1 

^ V 

x.n 


u, l,n ^u, l,n - l^n 


u,2,n '^u, 2,n " ^n 


W 


/3, = tan"l — 

T r 


x,n 


tan-l 

^x,n 


( 45 ) 

( 46 ) 

( 47 ) 

( 48 ) 

( 49 ) 


Critical velocity ratios are computed for the last stage, where temperatures are lowest 
and the velocity ratios are most severe: 




+ V" 


x,ex 


( 50 ) 
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In order to establish the flow annulus geometry near the turbine inlet, a flow analy- 
sis is made at the first-stage exit as follows: 


Ah’ =- 


gjx 


Ah’ 

T' = t: 

2, a in „ 


P2,a"Pin|l 


Ah’ 
a 


c T: 17’ , 
p in 'ay 


U 

AV = -^ 
u,a , 


V = -Hil AV 
^>2, a u,a 

u 


V 

„ = yT -Hii AV, ^ cot a, 

,^,a T u,a 1 

u 




a x,2,a 


T" _ >T<t 2, a 

2, a ■ ^2, a “77“ 
2gJc„ 


'T \r/(r-i) 

P2,a=P2,a(;jf‘‘’ 


■2, ay 


P2 a "" ^ 

^,a 


P2^ 

2, a 


(61) 

(62) 

(63) 

(64) 

(65) 

( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 
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( 71 ) 


A 


an, 2, a 


w 

^2,a^x, 2, a 



^an, 2,a 


7tD 


m,in 


1 , ^an,2,a 
m,in 


D 


t,2,a 



D 


h,2,a 



( 72 ) 


( 73 ) 


( 74 ) 


When a constant annulus is assumed for the first stage, the first-stage exit dimensions 
become the turbine inlet dimensions. 


Alternative Calculation Procedure 

The basic calculation procedure described in the previous section requires as inputs 
the inlet and exit mean diameters, stator exit angle, and shaft power. Alternatively , 
the hub or tip diameters could be specified as input and the mean diameters computed, 
the exit radius ratio could be specified as input and the stator exit angle computed, and 
the timbine pressime ratio could be specified as input and the shaft power computed. 

These alternative input options require iterative calculation procedures such as described 
in this section. 

With hub or tip diameters rather than mean diameters specified at the inlet and exit, 
it is necessary to assume initial values for the inlet and exit radius ratios. Initial values 
for inlet and exit mean diameters are then obtained as 
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(75) 


D = 
m 



when hub diameter is input and 


D 


m 



(76) 


when tip diameter is input. The computation then proceeds from equation (1) through 
equation (42) and the computed exit radius ratio is compared with the assumed value. If 
they are not the same (within a given tolerance), then the computed value of exit radius 
ratio is used to calculate a new value for exit mean diameter (from eq. (75) or (76)) and 
the computation procedure is repeated until convergence is obtained. Then, computa- 
tion proceeds through equation (72) and the computed inlet radius ratio is compared with 
the assumed value. If they are not the same, the computed value of inlet radius ratio is 
used to calculate a new value for inlet mean diameter, and the computation procedure is 
repeated from equation (1) . This entire procedure is repeated until both inlet and exit 
radius ratios in the same calculation pass converge to previous values. 

With exit radius ratio rather than stator exit angle specified, a value of stator exit 
angle is assumed for the evaluation of equation (16). The computation proceeds through 
eqviation (40) . Equation (42) is then used to compute the exit annulus area from the input 
value of radius ratio, and the exit axial velocity is then obtained from equation (41). The 
density used in equation (41), however, is not consistent with the exit area, and equa- 
tions (41), (33), (37), and (40) must be iterated until convergence is obtained. Then, the 
stator exit angle is computed as 


= tan-1 

V 

x.n 


(77) 


and compared with the assumed value. If they are not the same, the computed value of 
stator exit angle is used for the evaluation of equation (16), and the computation proce- 
dure just given is repeated until two consecutive values of stator exit angle are the same 
The remainder of the computation is then completed. 
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With the turbine inlet -total- to exit-static -pressure ratio rather than shaft power 
specified as input, an initial value of turbine static efficiency is assumed. Turbine work 

is then computed from 


Ah' 


nc t: 

' p in 



(78) 


instead of from equation (4) . The computation then proceeds through equation (26) in 
order to compute a static efficiency from equation (16). If the computed value is not the 
same as the assumed value, a new value of static efficiency is assumed, and the compu- 
tation is repeated until two successive values are the same. The remainder of the com- 
putation is then completed. 


DESCRIPTION OF INPUT AND OUTPUT 

This section presents a detailed description of the program input, normal output, 
and error messages. Included in the input and output sections are several example 
cases illustrating the use of the program and the various options. 


Input 

The program input, a sample of which is presented in table I, consists of a title 
card and the required physical data and option indicators in NAMELIST form. The title, 
which is printed as a heading on the output listing, can contain up to 77 characters lo- 
cated anywhere in columns 2 to 78 on the title card. A title card, even if it is left blank, 
must be the first card of the data package. Additional title cards can be used to identify 
different cases being run in the same data package. This is done by placing a title card 
in front of the data record for the particular case and using the option indicator ITIT as 
subsequently described. 

The physical data and option indicators are input in data records havir^ the 
NAMELIST name INPUT. The variables and indicators that compose INPUT and the 
proper units are as follows. These must be inputed for all cases except where otherwise 
indicated. Either the SI units or the U.S. customary units shown after them may be 
used. 
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PTIN 

TTIN 

MU 

R 

GAM 

DIN 

DEX 

RREX 


RPM 

POW 

W 

ALPHA 

KLOSS 

NMIN 

NMAX 

E 

PRTS 

lALPH 


IDIAM 


O O 

inlet total pressure, N/cm ; Ib/in. 

inlet total temperature , K; °R 

gas viscosity, (N)(sec)/m^; lb/(sec)(ft) 

gas constant, joules/(kg)(K) ; ft-lbf/(lbm)(^) 

heat capacity ratio 

inlet diameter - hub or mean or tip value as specified by the indicator 
IDIAM, cm; in. 

exit diameter - hub or mean or tip value as specified by the indicator 
IDIAM, cm; in. 

exit radius ratio; RREX may be omitted in the case where both IDIAM = 2 
and lALPH = 0; RREX is used as first trial value when lALPH = 0 and 
IDIAM = 1 or 3 

rotative speed, rad/sec; rpm 

shaft power - omit when IPR = 1, kW; hp 

mass flowrate, kg/sec; Ib/sec 

stator exit angle from axial direction; ALPHA is used as first trial value 
when lALPH = 1, deg 

turbine loss coefficient; a value in the range of 0.35 to 0.40 is usually 
applicable 

minimum number of stages for which the calculations are performed 

maximum number of stages for which the calculations are performed; 
results are obtained for all stage numbers between NMIN and NMAX 

squared ratio of stage-exit to stage-average axial velocities 

turbine inlet-total- to exit-static -pressure ratio; omit when IPR = 0 

indicates whether stator exit angle or turbine exit radius ratio is specified: 
lALPH = 0 - turbine is designed for specified ALPHA 
lALPH :: 1 - txirbine is designed for specified RREX 

indicates whether inputed diameters are hub, mean, or tip values: 

IDIAM = 1 - inputed diameters are hub values 
IDIAM = 2 - inputed diameters are mean values 
IDIAM = 3 - inputed diameters are tip values 
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IVD indicates type of velocity diagram used; 

IVD = 1 - symmetrical diagrams 
IVD = 2 - zero exit swirl diagrams 
IVD = 3 - impulse diagrams 

IVD = 4 - zero exit swirl diagrams if X ^ 0.5 and impulse diagrams if 
X < 0.5 

ITIT indicates use of title cards in addition to that required as first card of data 

package ; 

ITIT = 1 - title card precedes next data set; must be inputed for each addi- 
tional title card because ITIT is automatically restored to zero 
after each title card is read 

lEV indicates use of exit vanes: 

lEV = 0 - no exit vanes 

lEV = 1 - exit vanes are used to turn turbine exit flow to axial direction 

IPR indicates whether shaft power or pressure ratio is specified: 

IPR = 0 - shaft power is inputed 

IPR = 1 - turbine inlet-total- to exit-static -pressure ratio is inputed 

lU indicates type of units used for input and output: 

lU = 1 - SI units 
lU = 2 - U . S . customary units 

Each line of the input form shown in table I represents one data card. The first 
card is the mandatory title card, which can contain any desired message. The next 
three cards are the first data set, which contains all required inputs. This first case 
represents computation in accordance with the basic calculation procedure described 
previously. Data inputed for subsequent cases need only include those values that differ 
from previous case data. The fifth card is the second data set and represents the option 
where hub diameter is input. Also, the second case data specify that a title card, which 
is the sixth data card, precedes the third case data. Cards 7 to 10 represent four addi- 
tional cases illustrating use of different input options. The output corresponding to this 
sample input is described in the following section. 

Output 

The program output consists of title headings, the input variables, and computed 
results. This section presents normal output. Error message output is described in 
the next section. 
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Table II presents the output that corresponds to the sample input shown in table I . 
The top line of output is a program identification title that is automatically printed. The 
second line is the title card message. The next four lines are the input variables and 
their associated values for the first data set. The input variable names are spelled out. 
The units for the input variable values are as described in the Input section. The zeros 
printed under EXIT RADIUS RATIO and T-S PRESSURE RATIO indicate that these are 
computed for this case and not specified by the input. The fact that the input diameters 
are mean diameters is indicated by the MN in the variable name . Hub and tip diameters 
would be indicated by HB and TP, respectively. These four lines of output are printed 
for each new data set. 

The next two groups of eight lines each are the computation results for a one-stage 
tvirbine and a two-stage turbine, each satisfying the input requirements. Only one- and 
two-stage designs were specified by the input. The output parameters are spelled out 
and are self-explanatory . On the first line of each group are the number of stages, the 
stage speed-work parameter, and the diagram type, which is symmetrical for this first 
case. The remainder of the output includes exit and inlet tip and hub diameters in the 
first column, exit total and static temperatures and pressures in the second column, 
total and static efficiencies and velocity diagram angles in the third column, and first- 
and last-stage blade speeds and last-stage absolute velocity components in the last col- 
umn. The last line of each output group presents the last-stage absolute and relative 
critical velocity ratios. 

After the computations for each input case are completed, the input data for the next 
case are printed. The second case presented here is that where the specified diameters 
are hub values. The third input case in table I is preceded by an additional title card. 
This causes the next output to begin at the top of a new page with the pregram identifica- 
tion title and the title card message. The third to sixth cases are computed for two 
stages only. 

The third case is for a specified exit radius ratio rather than for a specified stator 
exit angle. This is indicated by a zero appeariig under STATOR EX ANG in the row of 
input variables. Exit vanes are included in this case, as indicated by WITH EXIT 
VANES printed after number of stages on the first line of result output. The fourth and 
fifth cases, as indicated by the top line of the fourth column of result output, are for 
zero exit swirl and impulse diagrams, respectively. Turbine pressure ratio rather than 
shaft power is specified in the sixth case. For this case, the inputed pressure ratio 
was chosen to be the same as the computed value obtained for the two-stage design of the 
second case. It is seen that, as should be expected, the program converges to identical 
solutions for both cases. An extra line of output consisting of shaft power and specific 
work is printed for the case where pressure ratio is specified as input. 
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Error Messages 


The program contains five output messages indicating the nonexistence of a solu- 
tion satisfying the specified input requirements. These messages are presented in this 
section, and their causes are discussed. 

(1) INSUFFICIENT ENERGY - This message is caused by the computed turbine exit 
total temperature being less than zero. It indicates that the turbine specific work re- 
quirement is greater than the energy available in the gas. Therefore, either the speci- 
fied shaft power must be decreased or the specified flow must be increased. 

(2) INSUFFIQENT IDEAL ENERGY - This message is caused by the computed ideal 
energy being more than that available from an infinite expansion of the gas. It indicates 
that the computed static efficiency is too low to yield a valid solution. Corrective action 
includes decreasing power, increasing flow, using more stages, or perhaps using a dif- 
ferent velocity diagram. 

(3) NEGATIVE TEX - This message is caused by the computed turbine exit static 
temperature being less than zero. It indicates a low value of turbine exit total tempera- 
ture and/or a high value of turbine exit velocity. Corrective action could be the same as 
for message (2) or decreasing a too high value of stage exit to average axial kinetic 

energy ratio. 

(4) INSUFFICIENT EXIT AREA - This message is caused by the computed exit area 
being larger than that available in the turbine. Such a situation can be remedied in many 
ways, including increasing exit diameter, decreasing stator exit angle, increasing inlet 
pressure, decreasing mass flow, and increasing stage exit to average axial kinetic 
energy ratio. 

(5) INSUFFICIENT INLET AREA - This message is caused by the computed inlet 
area being larger than that available in the turbine. Corrective measures are the same 
as for message (4). 


PROGRAM DESCRIPTION 

The computer program is called TURBAN. All computations are performed in one 
main program written in IBM 7090/7094 FORTRAN IV language. The program variables 
are defined in this section and the program listing is presented. 


A 

AA 


Program Variables 

-0 2 

factor in eq. (16), KRe * /cot 
factor in eq, (42), ^an,ex"^^^m,ex 
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AAl 


factor in eq. (72), A^^ ^ ^ 
turbine exit annulus area 


A EX 
AI 

AIL 

ALP 

ALPH 

ALPHA 

ALPHAl 

ALPHA2 

ALPHl 

ALPH2 

ASEX 

Al 

AIL 

A21 

B 

BETAl 

BETA2 

BETl 

BET2 

BL 

CCN 

CCP 

Cl 

CL 

CONV 

COT 

CP 

Cl 


2 

m,in 


stage loss parameter for intermediate stage and last stage havii^ no exit 
vanes 

stage loss parameter for last stage having exit vanes 

stator exit angle value printed in input section of output 

previous value of stator exit angle 

input value of stator exit angle 

output value of stator exit angle 

output value of stage exit angle 

stator exit angle 

stage exit angle 

speed of sound at turbine exit 

stage loss parameter for first stage or single stage having no exit vanes 

stage loss parameter for single stage having exit vanes 

first-stage-exit annulus area 

turbine exit loss parameter with no exit vanes 

output value of rotor inlet angle 

output value of rotor exit angle 

rotor inlet angle 

rotor exit angle 

turbine exit loss parameter with exit vanes 
dimensional constant 
dimensional constant 

blade row loss parameter for intermediate stage stator 

blade row loss parameter for exit vanes 

tolerance for radius ratio convergence 

cotangent of stator exit angle 

heat capacity 

blade row loss parameter for first-stage stator 
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D 

DELHT 

DELHI 

DEX 

DH 

DHEX 

DHID 

DHTID 

DHl 

DIN 

DM 

DMEX 

DMIN 

DN 

DT 

DTEX 

DTI 

DVUN 

DVUl 

DX 

D1 

E 

ES 

ESA 

ESI 

ESIL 

ESI 

ESIL 

ET 

ETI 


blade loss parameter for rotor 
turbine specific work 
first-stage specific work 
input value of turbine exit diameter 
data statement word HB for output use 
turbine exit hub diameter 

turbine ideal work based on inlet-total- to exit-static -pressure ratio 

turbine ideal work based on inlet-total- to exit -total -pres sure ratio 

first-stage-exit hub diameter 

input value of first-stage -exit diameter 

data statement word MN for output use 

turbine exit mean diameter 

first-stage-exit mean diameter 

turbine exit mean diameter 

data statement word TP for output use 

turbine exit tip diameter 

first-stage-exit tip diameter 

change in swirl velocity across last-stage rotor 

change in swirl velocity across first-stage rotor 

output word set equal to DH, DM, or DT as appropriate 

first-stage -exit mean diameter 

ratio of stage -exit to stage -average axial kinetic energies 

turbine static efficiency 

previous value of turbine static efficiency 

static efficiency of last stage with no exit vanes 

static efficiency of last stage with exit vanes 

static efficiency of single stage with no exit vanes 

static efficiency of single stage with exit vanes 

turbine total efficiency 

total efficiency of intermediate stage or last stage with no exit vanes 


ETIL 

ETl 

ETIL 

FLS 

FLSL 

G 

GAM 

lALPH 

IDIAM 

lEV 

IPR 

ITIT 

lU 

IVD 

J 

KLOSS 

LAM 

MU 

MXEX 

N 

NMAX 

NMLN 

NN 

PEX 

PI 

POW 

PRS 

PRT 

PRTS 

PTEX 


total efficiency of last stage with exit vanes 

total efficiency of first stage or single stage with no exit vanes 

total efficiency of single stage with exit vanes 

stator weighting factor 

exit vane weighting factor 

dimensional constant 

heat capacity ratio 

option indicator - see Input section 

option indicator - see Input section 

option indicator - see Input section 

option indicator - see Input section 

option indicator - see Input section 

option indicator - see Input section 

option indicator - see Input section 

dimensional constant 

turbine loss coefficient 

stage speed -work parameter 

gas viscosity 

turbine exit axial Mach number 
number of stages 
maximum number of stages 
minimum number of stages 
number of stages 
turbine exit static pressure 

TI 

shaft power 

computed value of inlet-total- to exit-static -pressure ratio 
computed value of inlet-total- to exit -total -pressure ratio 
input value of inlet-total- to exit-static -pressure ratio 
turbine exit total pressure 
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PTIN 

PT21 

P21 

Q 

Q2 

R 

RE 

REX 

REXl 

RHOEX 

RHOEXl 

RH021 

RPM 

RREX 

R1 

Rll 

ST 

STl 

ST2 

SUMUSQ 

TEX 

TITLE 

TTEX 

TTIN 

TTRN 

TTIN 

TT21 

T21 

U 

UISQ 


turbine inlet total pressure 
first-stage-exit total pressure 
first-stage -exit static presstire 

ratio of rotor inlet swirl velocity to change in swirl velocity 
ratio of rotor exit swirl velocity to change in swirl velocity 
gas constant 
Reynolds number 

computed value of turbine exit radius ratio 

previous value of turbine exit radius ratio 

turbine exit gas density 

previous value of turbine exit gas density 

first-stage -exit gas density 

rotative speed 

input value of turbine exit radius ratio 

first-stage-exit radius ratio 

previous value of first -stage -exit radius ratio 

output word set equal to STl or ST2 as appropriate 

blank data statement words for output use 

data statement words WITH EXIT VANES for output use 

sum of squares of stage blade speeds 

turbine exit static temperature 

input/output array for title card message 

turbine exit total temperature 

turbine inlet total temperature 

last-stage -rotor inlet and exit relative total temperature 

last-stage-stator exit absolute total temperature 

first-stage -exit absolute total temperature 

first-stage-exit static temperature 

rotor mean blade speed 

rotor mean blade speed squared 



VCRIN 

VCR2N 

VDl 

VD2 

VD3 

VD4 

VD5 

VD6 

VD7 

VD8 

VUIN 

VUll 

VU2N 

VU21 

VXN 

VXND 

VXl 

VXll 

VX2 

VX3 

VIN 

VIOVCR 

V2N 

V2NR 

V20VCR 

V21 

W 

WC 

WORN 

WUIN 


last -stage -stator -exit absolute critical velocity 
last -stage -exit absolute critical velocity 
data statement word SYMMET for output use 
data statement word RICAL for output use 
data statement word (blank) for output use 
data statement word ZERO E for output use 
data statement word XIT SW for output use 
data statement word HIRL for output use 
data statement word IMPULS for output use 
data statement word E for output use 
last-stage -rotor inlet swirl velocity 
first-stage-rotor inlet swirl velocity 
last -stage -rotor exit swirl velocity 
first-stage-exit swirl velocity 
turbine exit axial velocity 
last-stage average axial velocity 
first-stage exit axial velocity 

output word set equal to VDl, VD4, or VD7 as appropriate 

output word set equal to VD2, VD5, or VD8 as appropriate 

output word set equal to VD3 or VD6 as appropriate 

last-stage -stator -exit absolute velocity 

last-stage-stator -exit critical velocity ratio 

turbine exit absolute velocity 

last -stage -rotor -exit absolute velocity 

last -stage -rot or -exit absolute critical velocity ratio 

first-stage -exit absolute velocity 

mass flow rate 

factor in eqs. (56) to (59), [2ygR/(y + 1)]^^^ 
last-stage -rotor -inlet and exit relative critical velocity 
tangential component of last-stage -rotor -inlet relative velocity 
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WU2N 

WIN 

WIOWCR 

W2N 

W20WCR 

X 

XREX 

ZZ 


tangential component of last-stage-rotor -exit relative velocity 
last-stage -rotor -inlet relative velocity 
ist-st^e -rotor -inlet relative critical velocity ratio 
last-stage -rotor -exit relative velocity 
last-stage -rotor -exit relative critical velocity ratio 
function of y, r/(v - 1) 

exit radius ratio value printed in input section of output 
factor in eq. (35), 


Program Listing 


MBFTC TURBAN DECK 


C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


IVD=1 
IVD=2 
I V0=3 
I V0=^ 
IALPH=0 


THIS PROGRAM PERFORMS TURBINE GEOMETRY AND EFFICIENCY CALCS ON A 
MEAN SECTION BASIS ASSUMING SAME SHAPE DIAGRAMS FDR EACH STAGE 
lEXCEPT FIRST, WHICH HAS AXIAL INLET FLOW) AND IF 
SYMMETRICAL DIAGRAMS 
ZERO EXIT SwIRL DIAGRAMS 
IMPULSE DIAGRAMS 

ZtS FOR LAM. GE. 0.5 AND IMP FOR LAM.Lt.0.5 
EXIT RADIUS RATIO IS COMPUTED FOR INPUT VALUE OF ALPHA 

- INPUT RREX IS FIRST TRIAL VALUE IF IDIAM=1 OR 3 
IALPH=1 - ALPHA IS COMPUTED FOR INPUT VALUE OF EXIT RADIUS RATIO 

- INPUT ALPHA IS FIRST TRIAL VALUE 
ARE^INPUT AT INLET AND EXIT, AND BLADE SPEED VARIES LIN. 
INPUT DIAMETERS ARE HUB VALUES 

input diameters are mean values 
Input diameters are tip values 

NO EXIT VANES 

EXIT VANES TO TURN FLOW TO AXIAL DIRECTION 
POWER IS INPUT AND PRESSURE RATIO IS COMPUTED 
pressure RATIO(T-S) is INPUT AND POWER IS COMPUTED 
TITLE CARO PRtCEDES NEXT DATA SET 
SI UNITS ARE USED FOR INPUT AND OUTPUT 
U.S. CUSTOMARY UNITS ARE USED FOR INPUT AND OUTPUT 


DIAMETERS 
IDI AM= 1 
lOI AM=2 
IDI AM=3 
IEV=0 
IEV=1 
IPR=0 
IPR=1 
IT1T=1 
IU=1 
IU=2 


REAL LAM,NN, J,MXEX,KLOSS,MU 

n I M E NS ION U(99),TITLE(13),ST1I3),ST2(3),ST(3) 

NAMELI ST/INPUT/PTIN, TTIN,MU,R,GAM,DIN,DEX,RREX,RPM,POW,W,ALPHA, 

I KLOSS, I AL PH, NM I N, NMAX , ID1AM,E,IVD,ITIT»IEV,IPR,PRTS,IU 
DATA DH,DM,DT/2HHB,2HMN, EHTPZ ih ftHZFRD 

DATA VD1,VD4,VD3,VDA, VU5,VD6,V07,VD8/6HSYMMET,5HRICAL,1H ,6HZER0 

1,6HXIT SW,3HIRL,6HIMPULS,1HEZ 
DATA ST1,ST2/1H , IH ,iH ,6H WITH ,6HEXIT V»6HANES / 

PI= 3.1416 
98 WRI IE(6,100) 
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100 FORM at ( iH 1 f ^9X » 33HTUR8 I Nl VELOCITY DIAGRAM ANALYSIS) 
READ(5,99) TITLt 

99 F0RMAni3A6) 

WRITE (6,10i) TITLE 

101 FOR.MATIIH ,i3A6) 

ITI T=0 

1 REAU(5, INPUT) 

GO TO (95*96)iIU 

95 J=i, 

G = i. 

CCN=2. *PI 
CCP=1000. 

GO TO 97 

96 J=7/8. 

G=3^.17 

CCN=60. 

CCP-550. 

97 Rl=.9 
ES=.8 

ALP= ALPHA 

REX^RREX 

XREX=REX 

IF{ I ALPH.EU.l ) ALP=0.0 

IF (lALPH.EU.O) XREX=0.0 
IF(1PR,EQ.0) PRTS=0.0 
IF ( I PR.EQ. 1 ) P0W=0,0 
DO 1=1,3 
ST( I ) = ST1( I ) 

14 IF(IEV.EO.I) ST(I)=ST2(I) 

GO ro ( 13, 15, 17 ) , IDI AM 

13 OX = i)H 

DMEX=C l.+REX )/2.*UEX/REX 
DM1 ,! = { l.+Ri ) /2.*DIN/R1 
GO TO 19 

15 DX=UM 


DMEX=DEX 
DMIn=OIN 
GO TO 19 
17 DX = l)T 


DMEX=I l.+REX )/2.*UEX 
DMI.<(=( l.+Ri) /2.*0IN 
19 WRITEI6,109) 

109 FORMATdH ) 


no FORMATa26HO SHAFT MASS INLET INLET ROTATIVE INLET 
^ STATOR GAS HEAT GAS TURBINE AXIAL 

2 P,JWER FLOW TEMP PRESS SPEED ,A2,6H DIA ,A2r'' 

3S EX ANG CONST CAPAC VISCOSITY LOSS VcL SO ‘’pRESsJ 
‘+65X,5HRATItl,19X,5HRATI0, 14X, 19HC0EF RATIO RATIO/ 

6Ell.3.F7.3,,f8.3l ^ f '-2.f »-2 . F9.8 ,2f 8. 2 ,F8. 3 . 


ALPH1= alpha*. 017453 
X = GAM/ IGAM-i. ) 

CP= X*R/J 

IF(IPR.EO.I) Dh I D=CP*TTI N* ( I 
IF(lPR.EO.i) PEX=PTIN/PRTS 


-( i./PRTS)**! l./X) ) 
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1F( IP^.EQ.i ) GO TO 53 
OELHT=CCP*POW/W/ J 
TTtX= TTI N-OELHT/CP 
IF( ITEX.LT.0.0) GO TO 20 

53 N=NM1M 
CO.'JV = . 01 

37 IFlIU.EQ.l) 0N=UMEX/1C0. 

IF(IU.E0.2) DN=UMEX/12. 

IF(IU.EQ.l) 01=OMIN/1CO. 

IF(IU,E0.2) Ui=DMIN/l2. 

11 U( i )=PI*RPM*01/CCN 
U(.M)=PI*RPM*DN/CCN 
NN= FLOAT(N) 

IF(^.EQ.l) SUMUSO=U(N)*U(N) 

1F( i.EO.l) GO TO 3 
SUMOSQ= O.C 
DO • I=lfN 

U(II= (U(N)-Ol 1) )/(NN-l. )*(FLOAn I )-!.» + 0(1) 
UISO=J ( I ) **2 

2 SUMOS3=SUMU5Q+OISQ 

3 IF( IPR.EQ.O) GO TO 55 

54 DELHT=ES*DHi U 
TTbX=TTIN-DELHT/CP 
ESA=ES 

55 LAM^SUMUSO/G/J/DELHT 

GO TO (61, 6^., 63, 64), IVD 

61 0=(LAM+1. 

02=(LAM-1. )/2. 

FLS=2.-LAM 

VX11=VD1 

VX2=VD2 

VX3=VD3 

GO TO 65 

62 Q=1.0 
Q2= -.0 
ELS=1,0 
VXli=VD4 
VX2=VD5 
VXJ=VD6 
GO TO 65 

63 C=LAM+.5 
Q2=LAM-.5 
FLS-1. 0 

IF(LAM.LT..5) FLS=2. * ( l.-LAM ) 

VXll=tfD7 
VX2 =V08 
VX3=VD3 
GO TO 65 

64 1F(LAM.GE..5) GO TO 62 
GO TO 63 

65 DVU -i=J ( M) /LAM 
VUl ,= 0*DVON 
VU^;,4=32*DVUW 

4 COT= GDTAN(ALPHl) 

Cl=(l.+2.*COT**2)*Q*Q 
CI= Cl+(0-i. )**2 

0=2 . *C0T**2»Q*(J+ (U-L AM ) **2+ ( Q-LAM-1. »**2 
RE= W/MU/D1*2. 

IF(M.EO.l) KE=W/MU/DN*2. 



A=KLOSS/COT/RE**.2 
Al= A*( Cl+2. *0) 

AI=A*{ FLS*CI+2.*D) 

B= fc* CD T*»2»Q*U ■•-02**2 

ET1 = LAM/ ( LAM■^Al/2. ) 

ESi= LAM/(LAM■•■{Al■^B)/2. » 

ETl= LAM/(LAM■^AI/2, ) 

ESI= LAM/ ( LAM-M AI■^B) /2. ) 

IF (lEV.EQ.i) GO TO 3A 
1 F(m- 1) 5,t),6 
5 ES= ESI 
GO TO 7 

GO TO 7 

3A CL=2 .*COT**2*0*0+02**2 
FLSL=l. 

A1L=A1+A*CL*FLSL 
AIL=AI +A*CL*FLSL 
BL= E*C0T**2*Q*0 
£TlL=LAM/(LAM-fAlL/2. I 
ES1L = L AM/ ( LAM+ ( AlL-*-B L) / 2 . ) 
ETIL=LAM/ (LAM+AIL/2. ) 
ESR=LAM/(LAM + ( AIL+BL )/2.) 

3 5, j 5, 36 

35 ES=FS1L 
GO ro 7 

7 IF(IPR) 56,56,57 

56 ZZ=0ELHT/CP/TTIN/ES 

IF ( /Z. GE. 1.0 ) GO TO 21 
PEX=PTIN*( l.-ZZ )**X 
GO TO 56 

57 IF(ABS(ES-ESA).LT..0001) go to 53 
ES=(ES>ESA)/2. 

GO TO 5A 

58 VXND=0*COT*DVUN 
VXN=VXi^D*SORT ( E ) 

V2.''*= SORT { VU2N**2-*-VXN**2 ) 

V2NR = \/2N 

IF([EV.EQ,1) V2N=VXN 

TEX= TTEX-V2N**2/2./G/J/CP 

IFITEX.LE.O.O) go to 22 

IF( lU.EO.l) RHQtX=PEX/R* X 0000. /TEX 

I F ( I U. EO. 2 » RHOEX=PEX/R* .4 A. /TEX 

IF ( I AL PH. EO. i ) GO TO 8 

A£X= W/RHOtX/VXN 

AA= AEX/PI/UN**2 

IF ( AA. GE. 1 .0 ) GO TO 23 

REX .=RBX 

REX= ( l.-AA)/(l.+AA) 

IF( IDIAM.EQ.2) GO TO 9 

IF('J.EQ.l) CONV=.OO01 

IF ( ABS ( REX-RcXi ) . LT. CONV ) GO TO 9 

IF( 1DIAM.E0.3) GO TO 33 

DMEX=( 1 .+REx ) /2. *UEX/REX 

GO TO 37 
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3j DMtX=( l.-*-RLX)/2,*DEX 
GO ro 37 

8 AEX= PI*DN**2*( l.-REX )/{ L.+REX) 

31 VXN= rJ/AEX/RHQEX 

V2N=S0RT ( VU2N**2+VXN**2 ) 

V2NR = »/2N 

IF(IEV.EQ.l) V2N=VXN 
TEX^Tf EX-V2.'j**2/2./G/ J/CP 
RHO :Xl = RHOEX 

IF( lU. EO. 1 I RHOLX = PEX/R*i 0000. /TEX 

IF( (U.E0.2) RHDtX=PEX/R*:44./TEX 

I F ( A8S ( RHOEX-RHOEX 1 > .GT . .00 1*RHUEX » GO TO 31 

ALPH= ALPHl 

VXND = \^XN/SORT ( E ( 

ALPH1= ATAN2( VUIN.VXND ) 

IF( \BS( ALPHi-ALPH).GT..O 02) GO TO 4 
q ALPh2 = ATAN2( VU2N, VXND ) 

PTtX=PEX*(TTEX/TEX)**X 

DHT10=CP*TTIN*( 1 . - ( P T FX/ PT I N I ** ( 1 . /X ) ) 
ET=DELHT/DHT ID 
P0W=DELHT*K/CCP*J 
WUlN = y/UlN-U( N) 

WU2'J=)/U2N-U( N) 

HETI=ATAN2( WUIN.VXND ) 

BET ^=ATAN2(wU2N. VXND) 

W2N=SQRT(WU2N**2+VXN*«2 ) 

W1N = S0RT( WUiN**2 + VXND**2 ) 

V1N=S3RT( VU1N**2+VXND**2 ) 

TTR )=TTEX-( V2NR**2-W2N»*2 ) /2./G/ J/CP 
TTl J=TTRN-( W1N**2-V1N**2 )/2. /G/J/CP 
WC=SQRT(2.*GAM/(GAM+1. )*G*R) 
VCR2N=WC*S0RT(TTEX) 

WCRN=^C*SORT (TTRN) 

vcr;n=wc*sqrt{ttin) 

V10VCR=V1N/VCR1N 
WlOrtCR=WlN/WCRN 
W20WCR=W2N/WCKN 
V20VCR=V2NR/ VCR2N 
ASEX= SORT (GAM*G*R*TEX ) 

MXEX= VXN/ASEX 

DTEX= DMEX*2./ ( l.+REX ) 

DHEX= REX*OTEX 
IF( V.EO.l ) GO TO 51 
DELH1= U( 1 )**2/G/ J/LAM 
TT2..= TTIN-DELHi/CP 

PT2t= PTIN*l 1,-DELH1/CP/TTIN/ET1»**X 
DVUl = U(1)/LAM 
VX1=0*C0T*0VU1*SQRT( E ) 

VU21= Q2*0VU1 
VUll= 0*DVUi 

V21= SQRT(VU21**2+VX1**2) 

T21= TT21-V^l**2/2./G/J/CP 

P2i= PT21*(T21/TT21)**X 

IF ( lU. EO. 1 ) RH02 1=P21/R* 10000. /T21 

IF(IU.E0.2) RH021=P2 1/R* 144. /T2i 

A21=W/RH021/VX1 

AA1=A21/PI/D1/D1 

IF( AAl.GE.1.0) GO TO 24 

R11=R1 
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Rl= ll.-AAl )/ ( l.+AAl ) 

IF( lOI AM.E0.2 ) GO TO AQ 

IF( ABS ( Rl-Ri 1 ) .LT.CONV ) GO TO 41 

44 IF(IDIAM-2) 4;>,49,47 

45 DM1 l=( l.+Rl )/2.*DIN/R1 
DMtX = ( 1.+RL-X)/2.*DEX/REX 
GO TO 37 

47 DMig=( l.+Rl )/ 2 .*oiN 
DMEX=( l.+REX)/2.*DEX 
GO TO 37 

41 IFlCO'^V-.OOi ) 42»42f43 

42 GO rO 49 

43 CONV=.0001 
GO TO 44 

49 DT1=D'1IN*2./(1.+R1) 

DHl = R1*DT1 
GO TO 52 

51 R1=<EX 
DTl=OTEX 
DHl =DHEX 

52 ALPHA1=ALPH1/, 017453 
ALPMA2= ALPH2/.U17453 
BETA1=BET1/. 017453 
BETA2=BET2/. 017453 
PRT= PTIN/PTEX 

PRS= PTIN/PEX 

WRI TE(6,120)N,ST,LAM, VX1.,VX2,VX3.DTEX,TTEX,ALPHAI,U(1),DHEX,TEX, 
1ALPHA2,U(N) ,REX,PTEX,BETA1,VU1N,DT1,PEX,BETA2,VU2N,DH1,PRT,ET,VXN, 

2R1,PRS,ES,MXEX 

120 FOR'"AT ( 8H0STAGES = , 12 t 3 A6 ,2X, 27HSTAGE SPEED-WORK P ARAME TER = , F 5 , 3, 
126X,13HDIAGRAMS ARE , 3A6 

1 /20H EXIT TIP DIAMETER =tF6.2,4X, 18HEXI T TOTAL 

1 TEMP -,F7,2t5X,18HSTATUR EXIT ANGL E=, F6. 2 t 4X ,23HF I R ST STAGE MEAN 

2 SP-E0=,F8.i;/2GH EXIT HUB DIAMETER = t F6, 2» 4X t IBHEX I T STATIC TEMP = 
3,F7.2t5X,18HSTAGE EXIT ANGLE =, F6. 2 , 4X, 23HL AST STAGE MEAN SPEFU =, 
4F8.2/20H EXIT RADIUS RATIO = t F6 . 4, 4X, 18HEXI T TOTAL PRESS =,F7.^,5X 
5,18HROTOR INLET ANGL E = , F 6. 2 , 4X , 23HLAST STAGE INLET SW I RL = ,F 8. 2/2GH 
6 INLET TIP DIAMETER=,F6, ,4X, 18HEXIT STATIC PRE SS= , F 7, 2 , 5X , 1 8HROTO 
7R EXIT ANGLE = » F6 . 2 • 4X , 2 3HL A ST STAGE EXIT SWIRL =fF8.2/20H INLET H 
SUB OIAMETER=,F6.2,4Xt i8HT-T PRESS RATIO = t F7. 3 » 5X , 18HTOT AL EFFICI 
9ENCY = » F5 . 3f 5X t 23HLAST STAGE AXIAL VELOC=»F8.2/20H INLET RADIUS RA 
1TI0=,F6.4,4X, 18HT-S PRESS RATIO = , F 7. 3 t 5X, 1 8HSTAT IC EFFICIENCY=, 
2F5. 3,5X,23HEXIT AXIAL MACH NUMBER= , F 7. 4 ) 

IF( IPR.EQ.l. AND. lU.EQ.l) 0ELHT=DELHT/1000. 

IF(IPR.EO.l) WRITE(6,121) POW.DELHT 

121 F0RMAT(17H SHAFT POWER = . F9. 1 , 4X , 1 8HSPEC I FIC WORK =,F7.2) 

WRI TE( 6, 122 ) VlOVCRfWlQWCRf W20WCR, V20VCR 

122 F0RMAT(20H LAST STG ( V 1 / VCRl ) = , F6. 4, 4X, 18HL AST STG 1 Wl/ WCR 1 ) = , F 4 , 

15X,.8HLAST STG ( W2/WCR 2 ) = , F 6. 4, 4X , 23HL AST STG (V2/VCR2) =,F7.4) 

12 N=N+1 
CONV=.01 

IF( N.LE.NMAX ) GO TO 11 
IF( ITIT.EQ.i ) GO TO 98 
GO TO 1 

20 WRITE(6,130) 

130 FORMAT(1HO,5X,19H1NSUFFICIENT ENERGY) 

GO TO 1 

21 WRI TE( 6,140)N 
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lAO FORMAT ( IHOf 7HSTAGES= » 1 2 t 3X t 2 5HINSUFF IC I6NT IDEAL ENERGY) 

GO ro 12 

22 WRI fE( 5,150) N 

150 FORMAT ( lH0,7HSTAGfcS= , I 2, X , 12HNEGATI VE TEX) 

GO TO 12 

23 WRI TE(6,160)N 

160 FORMAT I IHO, 7HSTAGES= , I 2, 3X , 22HINSUFF ICI ENT EXIT AREA) 
Rl=.9 
REX=.7 
GO TO 12 

2A WRITE(6,173) N 

170 FORMAT ( IHO , /HST AGES= , I 2 , ^X , 2 3HI NSUFF IC I E NT INLET AREA) 

Rl = . 9 
REX=.7 
GO TO 12 
END 
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table n. - OUTPUT FORM WITH SAMPLE DATA 

TURBINE VELOCITY DIAGRAM ANALYSIS 
U.S. CUSTOMARY UNITS ARE USED. 
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TABLE n. - Concluded. OUTPUT FORM WITH SAMPLE DATA 

TURBINE VELOCITT DIASRAM ANALYSIS 
LED FOR HY ITIT IN PREVIOUS CASE DATA 
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Figure 1. - Effect of diagram type and speed-work parameter on velocity diagram shape. 
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